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Abstract. We discuss a numerical analysis employing the density of partition func- 
tion zeroes which permits effective distinction between phase transitions of first and 
second order, elucidates crossover between such phase transitions and gives a new 
way to measure their strengths in the form of latent heat and critical exponents. 
Apphcation to a number of models demonstrates the efficacy of the technique. 

A central theme of statistical physics is how best to distinguish between 
phase transitions of first and second order from simulational data for finite 
systems Numerical methods usually exploit the finite-size scaling (FSS) 
behaviour of thermodynamic quantities exhibiting rounded and shifted peaks 
whose shape depends on the order and the strength of the transition and 
which become singular in the thermodynamic limit at the transition point. An 
alternative strategy is the analysis of the FSS behaviour of partition function 
zeroes For field-driven phase transitions one is interested in the Lee- 

Yang zeroes in the plane of complex external magnetic field h and for 
temperature-driven transitions the Fisher zeroes in the complex temperature 
plane are relevant For d-dimensional systems below the upper critical 
dimension, the FSS behaviour of the j**^ partition function zero (for large j) 
is given by ||] 



Here, L is the system size, 77 is the anomalous dimension, t = T /Tc — 1 is 
the reduced temperature which is zero in the first formula, and h denotes 
the external field which is zero in the second formula. The integer index j 
increases with distance from the critical point. 

Despite early efforts , it has been considered prohibitively difficult if not 
impossible to extract the density of zeroes from finite lattice data [|j . This 
problem has resurfaced recently as zeroes-related techniques have become 
more widespread [Q. This provides the motivation for the present work in 
which we wish to suggest an approach suitable for density analyses. 

For finite L, the partition function may always be factorized as Zl{z) = 
A{z) (z — Zj{L)), where z stands generically for an appropriate function of 
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field or temperature and A(z) is a smooth non-vanishing function. Following 
Suzuki 1^ and Abe ||], we assume the zeroes, Zj{L), (or at least those close to 
the real axis and hence determining critical behaviour) are on a singular line 
for large enough L, impacting on to the real axis at an angle ip at the critical 
point z = Zc- The singular line is parameterised by z — Zc + r exp (iip). If we 
define the density of zeroes as (with Zj = Zc + rj exp (iif)) 

gL{r)^L-^'Y.^{r~r,{L)) , (2) 

j 

then the cumulative distribution is a a step function, 

GL{r)= f gL{s)ds=j/L^ if re(r„r,+i) . (3) 



It is natural to assume that at a zero, this distribution function is given by 
the average GUtj) = {2j - 1)/2L^. 

In the thermodynamic limit and for a phase transition of first order Lee 
and Yang Q showed that the density of zeroes has to be non-zero crossing 
the real axis, gooir) = 5oo(0) + a\r\'^ + . . .. For the cumulative distribution 
this implies the functional form 

G(r)=5oo(0)r + fe|rr+i + ... , (4) 

with the slope at the origin being related to the latent heat (or magnetization) 
via ^ goo(O) oc Ae. At second-order phase transitions, Abe and Suzuki [g| 
have shown that the necessary and sufficient condition for the specific heat 
to have the leading critical behaviour C ^ is gooir) = Ar^^" or 

G{r) cx r^-" . (5) 

For a = 0, as is the case in the d = 2 Ising model, it has been demonstrated 
[H that (^) leads to the correct logarithmic divergence in the specific heat. 

The preceding considerations show that a plot of Ghirj) — {2j — 1)/2L'^ 
against rj{L) should (i) go through the origin, (ii) display L- and j- collapse 
and (Hi) reveal the order and strength of the phase transition by its slope 
near the origin, parameterised generically as 

G(r) = air"' + as • (6) 

In order to test the efficiency of the density method we have examined six 



different models [11 for which various sets of partition- function zeroes exist in 
the hterature: 2D 10-state Potts, 3D 3-state Potts, 3D and 2D Ising, (3-f 1)D 
SU(3), 4D Abelian Surface Gauge. Here we shall summarize the results for 
the first three of these models. 



CD 



0.03 



0.02 



0.01 



0.00 
0.00 



Analysis of the Density of Partition Function Zeroes 



0.01 



0.02 

r 



0.03 0.04 



Fig. 1. Distribution of the L = 16-64 Fisher zeroes for the 2D 10-state Potts model. 
The symbols x , +, A, o, □ , and O correspond to j = 1 — 6, respectively. 



2D 10-state Potts model: The two-dimensional g-state Potts model is the 
classic testing ground for analytical and numerical studies of first-order [q > 
4) and second-order (q < 4) phase transitions. Apart from the transition 
point, /3o = ln(l -I- ^yq), also the critical exponents {q < 4), the latent heat 
{q > 4) and various other moments at /3o arc known exactly. 

Our analysis of the density of zeroes as tabulated in Refs. be gms 

with Fig. where the first six Fisher zeroes are plotted for L — 16-64. We 
observe excellent L- and j- collapse indicating that GL{rj) — {2j — 1)/2L'^ 
is the correct functional form. Fitting (||) to the j = 1-4 points gives 02 = 
1.10(1), 03 = 0.00004(1), strongly indicative of a first-order phase transition. 
Indeed, fixing 03 = 0, and fitting the two remaining parameters to the lowest 
four data points gives 02 = 1.008(6). Assuming thus 02 = 1, 03 = 0, and 
applying a single-parameter fit to the full data set we obtain g{Q) — ai — 
0.501(8). Further fits close to the origin yield the slopes and corresponding 
estimates for the latent heat Ae ~ 27rg(0) exp(— /3o) indicated in Table |l]. 



3D 3-state Potts model- 



It is generally accepted that this model exhibits 

This is there- 



a first-order phase transition, albeit a very weak one 1 14 



fore a typical model used to test new methods to discriminate between 
first- and second-order transitions. A list of the first few Fisher zeroes (with 
z = exp (-3/3/2)) for L = 10 - 36 can be found in Refs. ^Wi 



Table 1. Fits of the cumulative distribution G{r) to the A'' lowest Fisher zeroes 
(L = 16-64, j = 1-4) of the 2D 10-state Potts model. 
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9(0) 
Ae 


0.501(8) 
0.756(11) 


0.479(3) 
0.723(4) 


0.471(2) 
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Fig. 2. Distribution of the L — 18-36 Fisher zeroes for the 3D 3-state Potts model. 
The symbols o, □ , and O correspond to j = 1, 2, and 3, respectively. 



Our density analysis is presented in Fig. ^ A 3-paranieter fit to all data 
yields — 0.000005(2) and becomes even closer to zero as the fit is re- 
stricted closer to the origin. Clearly the slope is non-zero near the origin 
- the signal of a first-order phase transition. In fact, a 2-paraineter fit to 
the data corresponding to L = 22,24,30,36, j = 1 yields a2 = 1.06(2). 
Accepting that the plot is in fact linear near the origin, and fitting for the 
slope only gives 5(0) = ai = 0.0454(9). Using Po = 0.3670 ]ll,^3[|l| and 
Ae = 2TTg(0){3/2) exp(— 3/3o/2), we find that the corresponding latent heat is 
Ae = 0.247(5), comparing well with 0.2409(8) from ||l|,|ll and with 0.2421(5) 
from the more sophisticated analysis of ||l4| . 

3D Ising model: The first seven exact Fisher zeroes in the z = exp(— 4/3) plane 
for L — ^ are given in p5[ |, together with numerically determined zeroes for 
L = 5, j 1 - 4. We also use the zeroes in Rcfs. for L = 7, j = 1, 2; 

L = 6, 8, 10, 14, j = 1 - 3; i = 32, j = 1. 

Fitting the ansatz (||) to the full set of i = 4 — 32, j — 1 — 3 data 
shown in Fig. § indicates a second-order phase transition with a2 — 1.81(3), 
03 = —0.00001(1). Accepting 03 = and applying a 2-parameter fit to the 
six data points corresponding to L = 10 — 32, j = 1, gives 02 = 1.879(2) or 
a = 0.121(2), roughly compatible with the weighted "world average" ||lj of 
a = 0.10985(54). 

To summarize, we have shown that from the qualitative behaviour of the 
cumulative density of partition function zeroes we can distinguish between 
first- and second-order transitions while from the quantitative details we can 
extract the latent heat and the specific-heat exponent a, respectively. Our 
method meets with a high degree of success even in the borderline case of 
the 3D 3-state Potts model where with traditional methods the distinction 
between a first- and second-order phase transition is quite difficult. 
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Fig. 3. Distribution of the L = 4-32 Fisher zeroes for the 3D Ising model. The 
symbols x , +, □ , and O correspond to j = 1, 2, 3, and 4, respectively. 
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